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Abstract 

The plane wave decomposition method (PWDM) is one of the most popu- 
lar strategies for numerical solution of the quantum billiard problem. The 
method is based on the assumption that each eigenstate in a billiard can be 
approximated by a superposition of plane waves at a given energy. By the 
classical results on the theory of differential operators this can indeed be jus- 
tified for billiards in convex domains. On the contrary, in the present work 
we demonstrate that eigenstates of non-convex billiards, in general, cannot be 
approximated by any solution of the Helmholtz equation regular everywhere 
in (in particular, by linear combinations of a finite number of plane waves 
having the same energy). Prom this we infer that PWDM cannot be applied 
to billiards in non-convex domains. Purthermore, it follows from our results 
that unlike the properties of integrable billiards, where each eigenstate can 
be extended into the billiard exterior as a regular solution of the Helmholtz 
equation, the eigenstates of non-convex billiards, in general, do not admit 
such an extension. 



1 Introduction 

The quantum billiard problem in a domain Q C is defined (in units m=l) by 
the Helmholtz equation 

(_A _ k2)y,(a;) = 0, E = h\y2 (1) 
with Dirichlet boundary conditions 

^ix)\9n = 0. (2) 

The solutions En, ^Pn of these equations determine the energy spectrum and the 
set of eigenstates of Q. Studying the properties of {En, (pn) in quantum billiards 
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has became a prototype problem in "quantum chaos" . A simple form of eqs. ^ |21 
suggests a natural way to solve them. First, for a given energy E one looks for a 
set of solutions {tp^"'\k), n E IN} of the Helmholtz equation in the entire plane 
(without any boundary conditions). For example, {tp^^-^k)} can be chosen as a 
set of plane waves: {exp(ifc„a;), = k, fc„ G E^}, or as a set of radial waves: 
{ J„(kr) exp(m6')}. Then regarding {-ip^^^k)} as a basis one can search for solutions 
of eqs. m El using the ansatz 

y.(x)=5:a,^«(k,x). (3) 

As a result, solving eqs. [H El is reduced to the algebraic problem of finding the 
coefficients such that the linear combination vanishes whenever x G dQ. 

The above approach has been widely used both in analytical and numerical stud- 
ies of quantum billiards. In particular, it has been suggested by Berry in |Blj to use 
the expansion Q with a Gaussian amplitude distribution to represent eigenfunc- 
tions of quantum systems with fully chaotic dynamics. This idea has been applied 
in numerous works to calculate various quantities associated with eigenfunctions, 
e.g., autocorrelation functions |Blj . amplitude distributions |BSj . statistics of nodal 
domains |BGSj etc. The same strategy can be also used for a numerical solution of 
eqs. H El In this context it has been first introduced by Heller |Hej with the appli- 
cation to the Bunimovich stadium. Since that several modification of the method 
have been considered in |LR,Hj . |LHj and in |("LHj . Depending on the choice of the 
basis in the decomposition ^ one gets, in general, different numerical methods for 
solving eqs. dEl Here we will single out the basis of plane waves (PW), most often 
used in applications. For the sake of briefness we will refer to the corresponding 
numerical method as plane wave decomposition method (PWDM). 

As a matter of fact, the whole strategy described above is based on the as- 
sumption that the set {-ip^^^k)} furnishes an appropriate basis for the expansion of 
solutions of eqs. ^El In other words, one can use PWDM only if billiard eigenstates 
can be approximated by linear combinations of plane waves. That means 

W^n - ^t^^'lU^cn) ^0, as AT ^ OO (4) 
for some sequence of the states tp^^^ which are of the form 

N 

^[^]=5]a,e*^ hER", \h\=k. (5) 

i=l 

We will say that the plane wave approximation holds for a state if the limit P|) 
exists. 

Up to now it has been often assumed that the PWDM can be applied to billiards 
of arbitrary shape. From the results of Malgrange |Maj (see also [Ho]) on the theory 
of differential operators it is known that any solution of eq. Q regular in a convex 
open domain can be approximated by superpositions of plane waves with ki E C^. 
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Moreover, since each evanescent plane wave (Im/cj 7^ 0) can be approximated in a 
bounded domain by plane waves with real wavenumbers |B2j , one immediately gets: 

Proposition 1. Let Q C be a convex bounded domain, then any solution of eq. 
01 regular in Q can be approximated by plane waves. 

This shows that the eigenstates of a quantum billiard Q admit PW approximation 
inside any convex domain Qi C Q, see fig. la. Hence, PW approximation always 
holds for billiard eigenstates in a local sense. Furthermore, if is a convex domain 
one can choose Qi in such a way that dQi is arbitrary close to dQ. Consequently, 
as a simple corollary of Proposition 1 one gets: 

Corollary 1. Eigenstates of a convex billiard Q can be approximated by superposi- 
tions of plane waves. 

The question naturally arises whether the same property holds for eigenstates of 
non-convex billiards, and thus, whether the PWDM can be actually applied to the 
class of non-convex billiards. 

Note that there exists an important link between the PWDM and the problem 
of eigenstate extension in quantum billiards. Suppose y9„ is an eigenstate of Q which 
can be extended (as a regular solution of eq. from f2 to a convex domain Q2 ^ ^• 
Then it follows immediately by Proposition 1 that PW approximation holds for 
ifn- The example of a billiard where each eigenstate can be continued in a convex 
domain is shown in fig l.b. This is the "cake" billiard whose boundary consists of 
two concentric circle arcs connected by two segments of radii at an angle a < tt. In 
the polar coordinates x = (r, 6) the eigenstates of the "cake" billiard can be written 
explicitly as a sum of Bessel and Neumann functions: 

^^rr\^) = (C ^.Jk^r) + 6-njk(r)r)) sin {u^e - ^0)) , ^rn = ^. 

Since the singularity point of (p^^* (x) is always at the center O of the circle arcs it is 
possible to extend each eigenstate into a convex domain Q2, see fig l.b. Accordingly, 
any eigenstate of the "cake" billiard can be approximated by superpositions of PW. 

On the other hand, assume that for a billiard Q an eigenstate (fn can be ex- 
panded in a basis j-?/'*^")} (see eq. E}, where ip^^^^s are solutions of the Helmholtz 
equation regular in (e.g., plane waves). If furthermore, the corresponding sum 
(jni) converges everywhere in R^ it makes sense to consider ipn{x) both inside and 
outside fl. Such extension of (pn{x) into R^ provides simultaneously solutions for 
the interior Dirichlet problem (when x G fi) and for the exterior Dirichlet problem 
(when X E Q'^ = H'^/Q). Based on this observation a connection [spectral dual- 
ity) between the interior Dirichlet and the exterior scattering problems has been 
suggested by Doron and Smilansky in |DSj . The rigorous result has been estab- 
lished by Eckmann and Pillet |EPj . In most general form {weak spectral duality) it 



3 



could be stated as follows: En is an eigenvalue of the interior problem if and only 
if there exists an eigenvalue e~^^" of the exterior scattering matrix S{E) such that 
dn{E) —>■ 27T whenever E En- Moreover, if dn{En) = 27c {strong spectral dual- 
ity) then the corresponding interior eigenstate ipn could be extended into as 
functions. Therefore if strong form of spectral duality holds for some eigenenergy 
En then PW approximation holds for the corresponding eigenstate ipn- It has been 
explicitly shown that strong form of spectral duality holds for convex integrable 
bilhards |DiSj . However, as has been pointed out in [EPj . strong spectral duality 
cannot hold for billiards in general. 

Remark. It should be pointed out that the approximability by PW is much weaker 
property then strong spectral duality. As has been explained above, strong spectral 
duality implies PW approximation for the corresponding eigenstate. The opposite, 
however, is not true: PW approximation for an eigenstate does not imply, in general, 
strong spectral duality. In fact, in |B2t lEPj the examples of convex billiards (in 
this case the approximation by PW is possible) have been constructed where the 
eigenstates extension into the exterior domain as functions is not possible. 

2 Main results 

Let n be a simply connected bounded domain in R^ with a piecewise smooth bound- 
ary dQ. Two different billiard maps can be associated with Q. First, the standard 
billiard map \i/ corresponding to the motion of a pointlike particle in the interior 
domain. Second, the exterior map '^'^ which corresponds to the scattering off Q as 
an obstacle, see e.g., [Smj . In order to define the exterior map one can place i7 on a 
sphere of "infinite" radius. Then is a standard billiard map corresponding to 
the motion of a pointlike particle in the domain S'^/fl. It should be noted that there 
is an essential difference between convex and non-convex billiards. Whenever Q is 
a convex domain the interior map determines the same dynamics as the exterior 
map For any interior trajectory inside Q there is a dual trajectory in Q'^ which 
travels through the same set of points on the boundary dfl, see fig. 2a. We will refer 
to this property as interior- exterior duality. In particular, for convex billiards there 
is one to one correspondence between the interior and exterior periodic trajectories. 
For each periodic trajectory 7 its continuation 7^^ into the exterior domain will be 
the dual periodic trajectory of the exterior map. On the other hand, it is straight- 
forward to see that in non-convex billiards interior-exterior duality breaks down. 
Generally, in a non-convex billiard Q there exist interior periodic trajectories whose 
extension into the exterior domain intersects Q again, see fig. 2b. Let 7 be such a 
trajectory and let 7^^ be its extension in the exterior. Note that 7 U 7^^ is a union of 
straight lines in R^. Take / C 7 U 7^^ to be a line which intersects the boundary dfl 
at 2n, n > 1 points (for the sake of simplicity we will always assume that n = 2). 
Then the intersection 1] fl / is the union of two disconnected segments: 71 C 7 and 
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7i C 7^^. If 7i does not belong to any periodic trajectory in f2, we will refer to 
7 as single periodic trajectory (SPT). By definition any SPT has no dual periodic 
trajectory in the exterior domain. In what follows we call a non-convex billiard Q 
as generic if it contains at least one stable (elliptic) or unstable (hyperbolic) SPT. 
According to this terminology the "cake" billiard in fig. lb is non-generic, since all 
its periodic trajectories are of neutral (parabolic) type. 

We call a smooth function ijj{x) as a regular solution of the Helmholtz equation if 
it solves eq. ^ everywhere in R^. For a given energy E we will denote by Ai{E) the 
set of all regular solutions of eq. 1 and by A^pw(-E) C A4{E) the subset of functions 
which can be represented as linear combinations of finite number of plane waves 
with real wavenumbers ki, |/cjp = 2E/h'^. In particular, Ai{E) includes conver- 
gent superpositions of plane waves (also with complex wavenumbers i.e., evanescent 
modes) and radial waves with the energy E. In its crudest form the main result of 
the present paper can be formulated in the following way. Based on the breaking 
of interior-exterior duality we demonstrate that eigenstates of a generic non-convex 
billiard (in general) cannot be approximated by regular solutions of eq. H To illus- 
trate the main idea of our approach it is instructive to consider a non-convex billiard 
Q with an elliptic SPT 7. It is well known that a sequence of quasimodes (<^j,kj) 
associated with 7 can be constructed (see e.g., |PUlj . |Baj ). Each pair (<^„,k„) 
represents an approximate solution of eqs. Q El such that is localized along 7. 
Furthermore, in the absence of systematic degeneracies in the spectrum of Q the 
quasimodes ((^„, k„) approximate (in sense) a sequence of real solutions (v?n, k„) 
of eqs. m El For each such eigenstate (pn let us consider the corresponding Husimi 
function 

i/^„(z) = |(z|(^„)|2, z={q,p): qen, \p\ = hK, (6) 

where (z| denotes a coherent state localized at the point z of the phase space of 
Q. By the definition if<^^(z) is localized along 7 and exponentially small everywhere 
else. On the other hand, assume that ifn could be approximated by regular solutions 
of eq. m That means for any e > there is ipe ^ -^(-^n) such that | \ipn — iptW < ^, 
where || • || denotes the L^(f2) norm. Set g be a point at 71 and set p be directed 
along 7i. Then for z = {q,p) we have 

i7,„(z) =lim|(z|^,)r = lini|(z|e-^'^/Ve)r, (7) 

where e~**^/^ is the free evolution operator in R^. Furthermore, in the semiclassical 
limit the quantum evolution of coherent states is governed by the corresponding 
classical evolution 

e-*^/^|z) = e^*^/^|z(t)) + 0{h°^), z(t) = {q{t),p). (8) 

Plugging (jHl) into eq. Eland taking time t to be such that qit) = g' G 71 one gets 

H^A^)-H^,Xz') = 0{h^), z' = {q',p). (9) 
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This, however, contradicts the fact that the Husimi function H^^{z) should be ex- 
ponentially decaying outside 7. 

The above argument can be extended to the case of hyperbolic SPT 7 as follows. 
Contrary to the elliptic case it is not possible to construct quasimodes concentrated 
on hyperbolic periodic orbits. Instead, one can use a statistical approach in that 
case. By the results of Paul and Uribe |PUlj it is known that the average of the 
Husimi functions (jH)) 

{H,A^)) = ^ E (10) 

over the energy interval Vch = [E — ch, E + ch], c > depends in the semiclassical 
limit — *■ on whether z belongs to a periodic trajectory or not. On the other 
hand, as has been explained above, if each could be approximated by a regular 
solution of eq. [T]then each iJ^^(z) (and therefore the average (if^^(z))) would be 
(semi classically) invariant along 71 U 71. 

The preceding discussion provides an intuitive explanation why it is impossible to 
approximate eigenstates of a generic non-convex billiard by a superposition of plane 
waves. Speaking informally our argument says that contrary to the real eigenstates 
of non-convex billiard Q, any regular solution of eq. ^ always "preserves" interior- 
exterior duality. In what follows we consider the L^(f2) norm 

VnW = \\^r.-n, (11) 

for a solution (</)„, i?„) of eqs. [TJElin and an arbitrary ip G A4{En)- By the defi- 
nition rjniip) measures approximability of ipn by regular solutions of the Helmholtz 
equation. Recall that a state is approximable by PW if 

inf r]n{tp) = 0. 

Remark. Note that by Proposition 1 for any ip G Ai{En) and any e > one can 
always find G Mp^^En) such that |?7n(^) — Vn{4'e)\ < In particular this implies 

= , inf = inf Ui^). (12) 

In other words, an eigenstate can be approximated hj ip E J^{En) if and only if 
it can be approximated by PW. Therefore, in what follows one can always assume 
without lost of generality that belongs to Aipy{En) rather than to the set Ai{En)- 

By Corollary 1, 7]^'" = for any eigenstate of a convex billiard. On the contrary, 
in the body of the paper we show that for a generic non-convex billiard the average 
of 77^'" over an energy interval is bounded from below by a strictly positive constant: 
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Proposition 2. Let Vt he a non-convex billiard with at least one stable or unstable 
SPT and let ({pn,En), n = 1,2,...cxd denote the eigenstates and eigenenergies of 
the corresponding quantum billiard. For any set of approximating functions {ipi G 
Ai{Ei), i G IN} the average of rj^ = rin{ipn) over the energy interval Vch = [E — 
ch, E + ch], satisfies 

{rin)>C{h), where B=\imC{h)/h (13) 

is strictly positive and independent of ipi 's. Moreover, if Q contains a SPT 7 of 
elliptic type then (provided the spectrum of Q has no systematic degeneracies) there 
exists an infinite subsequence = {{ipj^, Ej^), m G IN} (of a positive density, i.e., 
liiriTv^oo ^'^■^"'1^'^^^ > such that for any G and any regular solution 

^ G M{En) 

r]^{ij)>C, + Oih'/'), (14) 
where is a strictly positive constant independent of ip and h. 

From p3|14j) one immediately obtains the corollary: 

Corollary 2. For a generic non-convex billiard Q there exists an infinite sub- 
sequence of eigenstates {fj^, n G IN} such that: 1) rfj")" > 0; 2) ipj^^ cannot be 
extended into the domain fl'^ (as a regular solution of eq. QJ). 

Obviously, this implies the following properties of a generic non-convex billiard: 

• In general, eigenstates of non-convex billiards do not admit approximation by 
PW and PWDM cannot be used in that case; 

• The spectral duality for a generic non-convex billiard holds only in the weak 
form. 

The paper is organized as follows. In the next section we collect several necessary 
facts about coherent states. In Sec. 4 the case of elliptic SPT's is considered. First, 
using the coherent states we construct a family of quasimodes ((^„, En) associated 
with such trajectories. Then, we show that the lower bound (|Tl|l holds for the 
eigenstates approximated by (^„. The case of hyperbolic SPT's is considered in 
Sec. 5. Here we use the results of Paul and Uribe to estimate the average (?7„) over 
an energy interval. Finally in Sec. 6 we discuss our results and consider possible 
generalizations. 

3 Coherent states 

Definition of coherent states. The coherent states have been introduced already 
in the beginning of quantum mechanics and have been used in many areas since 
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then. The basic idea is to built a complete set of vectors of Hilbert space localized 
in the phase space both in q and p directions at the scale \/fi. The standard example 
of such states in R'^ is given by the Gaussians: 



UI{X) = (det Ima)^ (/iVr)"^ ^i[{p,x-q) + l{x-q,a{x~q))] ^ ^ ^ ^^^^^^ > Q_ ^^5) 

In the present work we will consider a slightly more general class of coherent states. 
(For a more general definition of coherent states see e.g., jPUlj .) Let Pg(-) be a 
function in R*^ equal to one in a neighborhood of the point q and zero outside the 
sphere of radius e centered at q. A coherent state at z = (g,p) is the vector 

</)^(x)=p^,(x)<(x). (16) 
It is easy to see that the coherent states p6|) are semiclassicaly orthogonal: 

m\' = i + o{h), irM') = 0{h^) a z^z'. (17) 

The role of the cut-off Pq{x) is rather technical, it allows to define coherent states 
inside compact domains. To use the vectors (fT^ as coherent states inside a billiard 
domain Q one needs that 

supp[p^g(x)] C n. (18) 



Propagation of coherent states. An important property of coherent states is 
that their quantum evolution in the semiclassical limit is completely determined by 
the corresponding classical evolution. Let H = —h'^A/2 + v{x) be the operator of 
symbol Ti. = /2 + v{x) inducing the flow -.V -^V on the phase space V . Then, 
as it is well known, for any time t the propagation of the coherent state ip" localized 
at z G is given by 



e-*nA^. = ^i{sit)/n+m)^f) + 0{h^'^), (19) 



where S{t) = Jq (pq — H^p, q)) dt is the classical action along the path z(t) and /i(t) 
is the Maslov index. The parameters z{t) = \E'* ■ z, a{t) = ■ o in eq. ITOl are 
determined by the evolution of the initial data z, o under the flow \E'* : z ^ z(t) and 
its derivative 

D*' : . ^ .(() = (20) 

where d x d matrices a, 6, c, d are the components of D'^^ in a given coordinate 
system: 

I a h \ 



c d 



It is convenient to chose two of 2d coordinates in the phase space V to be along the 
flow and along the line orthogonal to the energy surface. Then the matrix a can be 
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decomposed into a = (t°©(t^, where the scalar part corresponds to the above two 
directions and {d — 1) x [d — 1) matrix corresponds to the orthogonal subspace. 
It is straightforward to see that in such a basis Z}\I^* acts separately on cx^ and a^. 
In particular, ■ a" is given by a linear transformation: 

D¥ ■ a° = f ■ (21) 

In the present paper we will use the above results for two types of two-dimensional 
flows: free evolution on under the Hamiltonian Hq {v{x) = 0) and the evolution 
induced by the billiard Hamiltonian Hq {v{x) = if x E fl and v{x) = oo otherwise). 
Let us consider in some detail the evolution of coherent states in billiards. Set fl 
be the billiard domain. We will denote by : V —>■ V the billiard flow, whose 
action is on the standard phase space V of Q. It should be pointed out that one 
can use the coherent states ()16|) for the point z = {q,p) eV only if q is sufficiently 
far away from the boundary dQ. Indeed, to satisfy the condition (fTH|) q has to be 
at the distance larger than e from the boundary. For the sake of simplicity, we will 
not consider a generalized class of coherent states defined in the whole domain Q, 
rather we will use the states (fT^ but only for the interior points of Q. For this 
purpose let us define the inner domain G Q which contains all the points q of Q 
such that the distance between q and dQ is larger then e: dist{q,dfl) > e, see fig. 
3. In what follows we will fix e to be a small compare to linear sizes of the billiard 
(but large compare to h^^"^) and consider the coherent states propagation under the 
condition that at the initial moment ti = and the final moment t2 = t the points 
z(0), z{t) belong to the domain fi^. Whenever this condition is fulfilled one can use 
the formula (fT^ . where the states 0^, (j)^^}^ are both of the form (fTBj) . Furthermore, 
if q{t) G VLe for all t G [^1,^2] (i-e., there is no colhsions with the boundary between 
the times ti and ^2) then the reminder term in (|T!?|l is of the order 0{h°°). 

Husimi functions. Let cpn be an eigenstate of H with the eigenenergy En- Given 
a coherent state 0^ one can construct the corresponding Husimi function: 

Hn{z) = z=(g,p); a = (a' , a'), -ta' = P > 0. (22) 

Based on the propagation formula (fTIH) the following average over Husimi functions 

00 771 771 

J2f{un)\Mr.)\' = j:di^~'''\ ^n = ^^, E=py2 (23) 

n 1=0 



has been calculated to the leading order by Paul and Uribe |PUlj . It turns out that 
the result depends on whether the classical trajectory through z is periodic or not. 
With the application to the Hamiltonian the results in |iPUlj read as follows. 
Let /(■) be the Fourier transform of /. If z is not periodic under the fiow then 




do= /(O). (24) 
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Alternatively, if z belongs to a periodic trajectory the additional terms (of the same 
order in h) arise. In particular, for a hyperbolic periodic trajectory 7 with the period 
the leading term in ()2H|1 is given by 

E /(^^t)— rT72777T ' (25) 




Z=— 00 



cosh^/^^;^^ 



where 

are the action, Maslov index and Laypunov exponent of 7. 

4 PW approximation for eigenstates of non-convex 
billiards (elliptic case) 

Let 7 be a periodic orbit in the billiard Vl and let T{E) be the "lift" of 7 to the phase 
space V at the energy E. This means T{E) is a set of the points z = {q,p) E V 
such that q E 'y, p"^ = ^E and the vector p is directed along 7. Obviously, for 
any z G r(i?), "if^ ■ z = z, where is the period of the trajectory. We will 
make use of the letter e to denote the restriction of 7, T{E) to the domain i.e., 
Y = {qe -fnile}, T%E) = {z = {q,p) E T{E) : q G l^J. Provided that 7 is elliptic 
a set of approximate solutions (quasimodes) (^„(x) of eqs. [H El associated with 7 can 
be constructed. The possibility of quasimode construction on elliptic periodic orbits 
is well known. In the following we will follow the approach developed in |PU2j . 
jPUlj (see also [Hcj, |Paj and the references there). 

Before we turn to the construction of the states ipn{x) in billiards let us recall a 
general definition for quasimodes. 

Definition. Let H he a Hilbert space and H be a self adjoint operator with the 
domain D{H). A pair ((^„,i^„) with (pn G D{H), \\ipn\\ = 1 and G R is called a 
quasimode with the discrepancy 5„, if 

{H - E„)ipn = Tn, with ||r„||=(5„. (26) 



By a general theory (see e.g., jLaj) the quasimodes (</?„, i?„) should be close to an 
exact solution {ipn, En) of the equation 

(H -E)(p = (27) 

in the following sense. If {ip, E) is a quasimode with the discrepancy 5 then there 
exists at least one eigenvalue of H in the interval 

Vs = [E-5,E + S\. (28) 
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Furthermore, let v be the distance between E and an eigenvalue Ei of H outside P^, 
then 

ll^-TT.^II < -, (29) 

where Hy denotes the spectral projection operator on the part of the spectrum {-En} 
inside the interval [E — i>, E + u). 

Remark. In general, the formula ()29p implies that any state cpn approximates a 
superposition of eigenstates (pn- In order to approximate individual eigenstates of 
H, 6n should be much less than the energy intervals: AEn = \En — -E'n+il, AEn-i = 
\En — En-i\- For two dimensional billiards {AEn) ~ so the approximation of ipn 
by (fn becomes semiclassically {h 0) meaningful only if the spectrum of Q has 
no systematic degeneracies and quasimodes with discrepancy 6 h", a > 2 can be 
constructed. For the quantum billiard problem a quasimode construction providing 
6 = 0{h°°) is known to exist |CPj and for the rest of this section we will assume 
that the billiard spectrum has no systematic degeneracies. 

4.1 Quasimode construction 

We will now schematically describe the construction of quasimodes concentrated on 
elliptic periodic orbits. The basic idea is to lunch a coherent state along the orbit and 
average over the time. As it can be shown, this procedure yields an approximately 
invariant state if the initial state is chosen in the right way, see e.g., |PUH I5cj. Let 
0^, z = {q,p) G r^(£') be a coherent state localized on the periodic orbit 7. We will 
associate with 7 the state 

\^hE)) = ^ ^''^"'-''''^^^ir.) dt, (30) 

where C is fixed by the normalization condition ||$r(£;)|| = 1 and is the period 
of the classical evolution along 7: z(T^) = z. The propagation formula (fT^ yields 

(E-H^,)<|.?(^) = r„ 

Cr^ = z^(e'(^-/^+'^-Vz^^"^ -C) +0(/i^/'), (31) 

where S^, fi^ are the classical action and Maslov index after one period. Therefore, 
Cr^ = 0{h^^'^) provided that the following conditions are satisfied: 

Condition 1: cr(T^) = a; Condition 2: S^/h + /i^ = 27rn for some integer n. 

For each n let £n^ Cn = (c"°, cr^) denote solutions of Conditions 1, 2. It is possible 
to show (see e.g., |PUlj ) that the first condition can be satisfied if and only if 
an = and 7 is an elliptic periodic orbit. The second condition impose the Bohr- 
Sommerfeld quantization on the quasienergy Sn- When both conditions are satisfied 
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the corresponding pair {Sn,^Yle„)) Provides the quasimode with the discrepancy 

= o{h^^^)/c. 

Remark. It should be noted that a much wider class of quasimodes concentrated 
on 7 can be constructed by this method if one uses in (jHUj) coherent states with 
transverse excitations |PUH IPaj . For simplicity of exposition, we restrict our con- 
sideration only to the quasimodes without transverse excitations, whose leading 
order is determined by eq. EDI 

To construct quasimodes with discrepancies of higher order in h one has to con- 
sider the time evolution of coherent states of a more general type. This leads to 
transport equations whose solvability pose additional conditions on the quasiener- 
gies, see [Hcj- From the results of Cardoso and Popov J3P\ the prossibillity to 
construct quasimodes {En, 'fn) in billiards having discrepancy 6^ = 0{h^) is known 
to exist. Let {s,y) be a coordinate system in a neighborhood of 7 such that s is 
a coordinate along the trajectory and y is a coordinate in the orthogonal direc- 
tion. Using these coordinates the leading order of (i?„, (pn) can be written as follows 

[EillSEl: 

En = £n + 0{h^), = e'^^^^/'^uix) + 0{h), (32) 

where 

vis, y) = vo{s)y^ + 0{y^), u{s, y) = uo{s) + 0{y^) 
and the parameters fo(s), uo{s) are determined by Conditions 1, 2: 

$-^,^)(x) = e»«(^)^'/^«o(s), x = {s,y). (33) 

As has been explained before, in the absence of systematic degeneracies in the 
billiard spectrum one can expect that, in general, a state cpn approximates an in- 
dividual eigenstate of the billiard Q. In what follows we will denote by the set 
of quasimodes for which approximates some eigenstate (pn (rather than a linear 
combination of V9„'s) and by the set of true solutions of eqs. HI El corresponding 
to S^. Then by eq. 1221 for each ((^j, Ei) G and {ipi, Ei) G we have 

C] = \\^i-^,\\=0{h'^), \E,-E,\ = 0{h'^). (34) 

4.2 A lower bound for the approximation of eigenstates 

The quasimode construction described in the previous section is quit general and 
can be applied to an arbitrary elliptic periodic trajectory. In the present section we 
will consider eigenstates of the billiard VL from the subset S^, where 7 is an elliptic 
SPT. We show that for (^9^, En) G and any regular solution ip G J^{En) of eq. Ill 
in the norm 

Vnii^) = \Wn - i^W (35) 
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is bounded from below by 

VnW>C, + Cl, + 0{n'/'), (36) 

where Cj is a positive constant determined only by geometrical parameters of the 
periodic orbit. Since = 0{h°°), this implies the inequahty (fT^ holds for any 

{ifn, En) E S^. 

Let 7 be an eUiptic SPT and let 71, 71 be as defined in Sec. 2, see fig. 3. 
Now fix the parameter e to be sufficiently small such that 7^ = 71 fl 7^ 0, 
^f = ^inn/^ 0. We will denote by the capital letters Ti{E), Ti{E) (resp. (E), 
Tl{E)) the corresponding "lifts" of 71, 71 (resp. 7f, 7f) into the phase space V at 
the energy shell E. Recall that the main idea behind the quasimode construction 
flHUjl is to use coherent states propagating along a periodic orbit. By analogy, one 
can construct states localized on 71 and 71. Let z(0) = z G Ti{E). Consider 
the classical evolution (both for positive and negative time) of z under the free 
flow \1'q : z — > z(t) = {q(t),p(t)) in R^. Obviously, as time evolves, the point q(t) 
successively crosses the boundary of Qs at the sequence of points (?2, ^1, ^2, see fig. 
3. We will denote by ^1,^21 ^1,^2 the corresponding time moments: qi = q(ti),q2 = 
9(^2)5 Qi = '?(^i), ^2 = lih)- Then the states localized along 71 and 71 are given by 

m,iE)) = dt, (37) 

\n,iE)) = iy^'^''-'''^^'m dt. (38) 

Note, that under the free evolution e"**'^"/^ the support of 0^ is not preserved inside 
fl, and therefore the supports of , do not belong to the billiard domain. 
However, one can slightly modify the definition of the states $f ^ , to make them 
admissible as billiard states in Q. Let z = zi, a = ai be as before and set r be such 
that under the classical evolution \E'q : zi — »• z(t) the point z(r) = Z2 belongs to . 
Set 0z2^(a;) = e~*'^'^''''^0^(x) + 0{h°°) be the coherent state in fi, whose parameters 
are given by: ((T2, 7.2) = {D'^q ■ cti, \E'5 ■ zi). Then the states 

m.(E)) = £e^*(^-^")/l0-)rft, (39) 

m.iE)) = r"^^*(^-"")/^i0-)rft, (40) 

Jt2—T 

have their supports in Q and satisfy 

m,iE)) = m.iE)) + oih^), = \n,iE)) + oih^. (41) 

To get the lower bound (|36|) we are going first to construct a state $ with the 
property 

(^1$) = + O(/i~), (42) 
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for any ip G M.{E'). Let us show how such a state can be constructed using $p^, 
Set (■|-)r2, ("I-) be the scalar products in L^(R^) and L^(f2) respectively. From 
the definitions ()H7|H8j) one has 

(V'|$?,(^))r2 = r e^*(^-^')/^(^|0n dt = C^{u){m), (43) 



t2 



where 

^/ X + t2)w\ 2sin(wT,,/2) ^ , , 
Ci(o;)=exp(^^y^ ^7^' ^71 = 1^1-^21, (44) 



and oj = {E — E')/h. Analogously: 

mhE))n^ = C2H{m), (45) 

with 

^/ X + t2)a;\ 2sin(cuT;.,/2) ^ _ , 

C2(cu) =exp ( -^f^, ^71 = 1^1-^21. (46) 

Furthermore, let us introduce the states 

Then it follows immediately from eqs. |33llSthat the state $ = ^'^{E,E'), 

E')) = E')) - miE, E')) (48) 

satisfies orthogonality condition (jl^ . 

Let i?„) G 5^ be a solution of eqs. [H |21 and let -E„) G 5^ be the 
corresponding quasimode, whose leading order parameters £n, cr^ = (cr°, cr^) are 
determined by Conditions 1, 2, see eqs. ESI Now fix the energy parameters in 
eq. EHl by = £n, E' = En and put a = o"„, where cxn = {ij3, a^) and j3 is an 
arbitrary real positive number. We will make use of the state 

in order to get a lower bound on rjn- For any G A^(i?„) we have 

W^n-^PW ll^nll > |(^n-^|$n)| = | (<^„ | $„) | + O (/l°^) . (49) 

Using the triangle inequality 

W^Pn -iPn\\ + llv^n " > 1 1 <^n " V^n + V^n " 1 1 = jjf^n " ^|| (50) 

one gets immediately from ()49j) 

r/.(^) = llv.. - ^11 > - + 0(;.-). (51) 

1 1 1 1 
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It remains to estimate the scalar product |(<y9„|$„)| and the norm of the state 
First, consider the norm ||$n||- Since 71 fl 71 = one has from the definition of 



with Un = {En — Sn)/h. The calculations of the scalar products performed in 
Appendix give: 



1/2 



and for the leading order of Ci(ci;„), C2{uJn) one has from eqs. 

\C2{uJn)\ = T^, + 0{n), \Ci{uJn)\ = T^, + 0{h). (54) 
Combining (I^Hj) and together one finally gets 



In the same way for the scalar product ((y9„|$„) we have by 



The estimation (j36|) follows now immediately after inserting eqs. EHlEElinto ()5H1 . 
The resulting constant C^, which determines the lower bound on rin in the semiclas- 
sical limit reads as 



Cry 



T- T 



^71 '^71 



+ 0{e), (57) 



where , , are the lengths of 71 , 71 and 7 respectively. 

5 PW approximation for eigenstates of non-convex 
billiards (hyperbolic case) 

In the present section we consider the case of a hyperbolic SPT 7. Let as before 
{iPri{.x)} be the set of eigenf unctions in VL approximated by regular solutions {■?/'„ (x)} 
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of eq. ^ For an arbitrary set of ipni^) G A4{En), n = 1, 2, ... oo we will estimate the 
average of 

Vn = Vnii^n) = WVn - 1pn\\ (58) 

over an energy interval. Our objective is to show that independently of the choice 
of ipnS, in the limit h the average (rin) is bounded from below by a strictly 
positive constant. 

Let ^i{E,E'), $2(-E',-E')) ^"{E^E') be as in the previous section with the pa- 
rameter a of the form a = {if3, a^), f3 > 0. For each integer n we will consider the 
states 

|$„,i) = \^UE, K)), |$„,2) = miE, En)) (59) 
and their difference 

\^n) = |$n,l) - |$n,2) = {E , K)), (60) 

which is orthogonal to any ip G M.{En) up to the term 0{h°°) (see eq. In 
addition, it will be also useful to introduce the state 

1$;) = |<l>„,i) + |<l>„,2). (61) 

Note that $^ is orthogonal to in the semiclassical limit. 

Similarly to the case of elliptic SPT's, one can make use of the state to get 
a lower bound on 77^: 



Vn > rr^—. = ||~ II h 0[h ). (62) 















1$ 1 

1 n \ 






n \ 





In order to estimate the right side of this inequality let us consider the following 
difference 

Vn = \ {^n,lM\'' - \ {^n,2M\^- (63) 

Using the states $„, one can rewrite Vn as 



Vn = Re[{<!>r.\^n){^'JiPny). (64) 
Hence, the following inequality follows immediately 

\Vn\ < |($n|^n)| < ||<|| |($n|^n)|. (65) 

Finally, since ||$„|| - ||$;|| = 0{h°°), we get by ^ and ^ 



+ 0{h°°). (66) 



We will now use this inequality to get a lower bound for the sum of r]n over the 
energy interval Vch = [E — ch,E + ch], where c is a positive constant. One has 
straightforwardly from 
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J2 Vn> 



E 



\{^n,l\fn) 



E 



\{^n,2\fn) 



(67) 



Furthermore, the definition of the states $n.i! imphes 

\{<^n,iM\' = miMl', zi e Tl; rnrMl' = KCl^n)!', G Tf, (68) 

where (zi,cri) = (z, cr) and (z2,o"2) = (z(r),cr(r)) are related by the free classical 
evolution as in the previous section. As a result, the inequality (jFTTj) reads as 



with Un = {E — En)/h and 

/(t^n 



E/WK^n|0^,OP-E/WK^nlC)l' 



+ 0(/i°°), (69) 



ifcU„G[-C,c] 

otherwise. 



The elementary calculations (see Appendix) provide the leading order of the function 

f{uJn), UJn e [-C, C]: 



+ o(r 



2|p| 



[nhpy^ \sm\ur,T^j2) sin2(^„T^^/2) ^ 



+ 0(^°). (70) 



Now we can apply to (jH^jl the results of Paul and Uribe (see Sec. 3). Taking 
into account that zi G F while Z2 does not belong to any periodic trajectory, we get 
by eqs. I^UHlthe following estimation for the average of 



1^0 



jil{Sj /h + fj.-y) 



where F[-) is the Fourier transform of the function 



cosh^/2(^^^ 



+ 0(/i=^/'), (71) 



F{x) 



8 \ 2 



sin^ {xTj-^ /2) ~'~ sii?{xT^^/2) 





if X G [— c, c] 
otherwise 



and i^Vcu is the number of eigenstates in the interval Vcn whose leading order for a 
billiard of area A is given by the Weyl formula: 



Consequently, if 



Y 



^ pil{S~f/h+fj.^) 

J2f{it^) — ^0 

i^Q cosh / (/A^' 



(72) 
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one has from ()71|) 

(r/„) >Bh + 0{h^^^), B = 27tY/cA>0. (73) 

If moreover one assumes that T^^^c, T-y^c << 1, the function F{x) takes a simple 
form: 

[ otherwise 

and the constant B can be written exphcitly: 

A \ T;yj + T^j J 

Note, that the lower bound (fTSj) has been obtained using only one SPT. In the 
case of hyperbolic dynamics, however, the periodic orbits (and, in particular, SPT's) 
proliferate exponentially. Therefore, one can improve the estimation (ffH|l making 
use of a state which is concentrated on a set of SPT's {7} and satisfies eq. |121 
A simple way to construct such a state is to define it as the superposition: 

^r = E'5n(7), (75) 
{7} 

where $„(7) stands for the state (plj) associated with a SPT 7. 

Finally, let us mention that the statistical estimation (ffH|l can be straightfor- 
wardly generalized to the case of elliptic SPT's. In that case one should use the 
analogs of eqs. |211 123 (which are known to exist jPUlj ) for stable periodic trajecto- 
ries. 



E 



sin (/cT^ 



COS. 



(74) 



6 Discussion and conclusions 

Speaking informally. Proposition 2 implies that there is no on-shell basis of regular 
solutions of the Helmholtz equation which can be used to approximate all eigenstates 
of a generic non-convex billiard. That means any linear combination of plane waves, 
radial waves etc., with the same energy fails to approximate real eigenstates of 
non-convex billiards. In fact, a stronger result can be shown. Let f2 be a generic 
non-convex billiard and let Q' be a domain (not necessarily convex) which properly 
contains Q: Q' D Q, dQ' fl dQ = 0. Denote by Ai^'{E) the set of all solutions of 
eq. HI regular in Q' (note, that M.qi{E) D A4{E)). Let us argue that the eigenstates 
of Q cannot be approximated, in general, by states belonging to AiQ'{E). Let 

7 be a SPT and let /,7i,7i be as defined before. Furthermore, assume that the 
segment of the line / between 71 and 71 is entirely in Q', see fig. 4. (It seems 
to be a natural assumption that in a generic case one can always fined such a 
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SPT, provided Q' properly contains Q). Then take f2o C to be a convex domain 
satisfying: Qq Hji ^ 0, Qq H •ji 7^ 0. Now, suppose an eigenstate ipn of Q can be 
approximated by states ip'{x) from A^n'(-^n)- According to Proposition 1 ip'^x) can 
be approximated in by regular solutions of eq. ^ and thus for any e > there 
exists ipt, G A4{En) such that — ipe{x)\\L2(^Q^^ < e. Therefore, applying the 

same arguments as in Sec. 2 we get 



where zi = G Ti{En), gi G 71 n Qq and Z2 = (g2,p) e ^i{En), g2 e 7i H fio- 

However, as has been pointed out before, this cannot be true for each n since Z2 ^ F. 

The two properties of generic non-convex billiards follows immediately from the 
above analysis. First, it is not possible to approximate eigenstates of a generic non- 
convex billiard f2 also if one includes in the basis {^^^\k)} singular solutions of eq. 
^ e.g., the Hankel functions 



with a finite number of singularity points Xj. Second, there exists an infinite se- 
quence of eigenstates which do not admit extension into any large domain Q' prop- 
erly containing Q. That means the continuation of the interior eigenstates of a 
generic non-convex billiard into the exterior domain should be (in general) impos- 
sible because of singularities which occur arbitrary close to the billiard's boundary. 
It remains as an open problem what is the exact nature of such singularities. (For 
example, whether one can, in principal, extend eigenstates beyond the boundary 
of a generic non-convex billiard.) It should be also mentioned that the problem 
of the eigenstates extension in convex billiards is beyond the scope of the present 
paper. It would become a natural question to inquire about the relation between 
the billiard shape and the type of singularities arising for the extended eigenstates. 
In particular, it would be interesting to know whether the strong form of spectral 
duality (when it is possible to extend eigenfunctions in as regular solutions of 
the Helmholtz equation) holds exclusively for integrable billiards. 

Further, let us stress an important difference between the cases of elliptic and 
hyperbolic dynamics. The counting function A/'*(k) = #{k„ < k} for quasimodes 
[ifn, k„) which can be constructed on an elliptic periodic trajectory is known to be 
of the same asymptotic form A/'*(k) = ak^ -|- 0(k), a > as the counting function 
A/" = Ak^/An + 0(k) for the real spectrum {k„}, see |(yPj . Therefore, in a generic 
case, if an elhptic SPT 7 exists the subsequence {'^j„,n G of billiard eigenstates 
approximated by the quasimodes concentrated on 7 should be of the positive density: 



Since for each cpj^ the estimation (fT^ holds, that means there exists a subsequence 
of eigenstates with a positive density which do not admit approximation by plane 



H^„iz^) - H^^iz,) = lim|(zi|z^,)|2 -lim|(z2|^.)r = Oih^^) 



n G IN}, 
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waves. In the case of hyperbolic dynamics the statistical lower bound ()13|) implies, 
in fact, only a weaker result. It says that an infinite sequence (possibly of zero 
density) of such states exists. However, if one assumes that all eigenstates of fully 
chaotic billiards have "uniform properties" the inequality ()13|) suggests a natural 
conjecture: 

Conjecture. For a non-convex billiard with fully chaotic dynamics the set of states 
which can be approximated by PW is of density zero. 

Note, that it is impossible to exclude the possibility of existence of "exceptional" 
eigenstates (the eigenstates which can be approximated by PW) in non-convex bil- 
liards. Indeed, one can take a finite superposition of plane waves tp^^^ and set a 
(non-convex) nodal domain of ^/'t^' to be the billiard's boundary. Then tjj^^^ itself is 
the eigenstate of this billiard which can be approximated by PW. 

Finally, the study of the present paper is restricted to the two-dimensional sim- 
ply connected domains with Dirichlet boundary conditions. However, it is easy 
to see that presented results allow several rather straightforward generalizations. 
First, higher dimensional billiards and different types of boundary conditions can 
be treated in the same way. Second, billiards in multiply connected domains (fig. 5) 
have the same properties as non-convex billiards. Consequently, all the results ob- 
tained for non-convex billiards hold for multiply connected billiards as well. Third, 
we conjecture that our results can be generalized to the billiards on non-compact 
manifolds with non-trivial metrics (also in the presence of a potential) e.g., billiards 
on the hyperbolic plane. In such a case, one needs to adjust the notion of domain's 
"convexity" to the corresponding classical dynamics. In other words, a domain 
should be defined as "convex" if the interior-exterior duality holds and defined as 
"non-convex" if it breaks dawn. 
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Appendix 

Proposition 3. Let $p, $p be the states: 




(76) 
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localized along the path T = T{E), T{E) = {^*-z = {q{t),p{t)), t E [0,T], E = p"^ /2} 
with cr" = iPi, ct° = i/32; l3i,P2 > and = a^. Then 

($-|$^) = ($^|$^)i/2($^|$^)i/2 + o{h). (78) 



Proof. The inner product 

C1C2 -/o 



($-|$^) = -J— r r(0^|e'(^-^o)(*^"*^)/^|0^)dtidt2 (79) 



can be written as 



J\t -\t\)H{t)dt, (80) 



1 



2C1C2 
where 

i7(t) = (0^|e'(^-^°)*/''|0^). (81) 
By the propogation formula (fTIHl we get for (jHT|) 

i/(t) = e'(^W+^)/"+*'^W(0^|0^g) + O(^) 

^ ((fS^-w) -(-TOO -'^">- 

After inserting this expresion into eq. IHOland applying the stationary phase approx- 
imation to the integral one gets ()77I78|) . Finaly, let us note that eq. ITSl remains true 
also when /?i or P2 equals zero. 
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